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Abstract 

This work studies a simplified model of the gravitational collapse of a protostar 
based on the assumption that the mean velocity of the gas is proportional to the lo- 
cal acceleration. This assumption leads to a nonlocal generalization of the restricted 
Patlak-Keller-Segel model. We discuss some of its fundamental properties and illus- 
trate them by means of numerical simulations. In particular, we are interested in 
the conditions for the gravitational stability of an infinite gas cloud, represented by 
T d . Local well-posedness in Sobolev spaces is proven for non-local pressure forces. 
Moreover, we show the smoothing effect of our equation, as well as a Beale-Kato- 
Majda-type criterion in terms of ||p||i<». We also show that the problem is ill-posed in 
Sobolev spaces when it is considered backward in time. Finally, we prove that, in the 
critical case (one conservative and one dissipative derivative), ||p||z,°° (t) is uniformly 
bounded in terms of the initial data for sufficiently large pressure forces. 
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1 Introduction 

Gravitational collapse - the infall of a body under its own gravity - is of paramount im- 
portance in Astrophysics. Stars and galaxies grow from initially small preturbations of 
a uniform density medium, which is an unstable equilibrium solution of the Newtonian 
equations of motion. Upon the action of gravity, overdense regions attract the neigh- 
bouring material, accreting all mass within a certain radius of influence into a singular 
point. Some of the gravitational energy, though, transforms into heat and/or random mo- 
tions that oppose collapse. In most real- life situations (black holes being the exception), 
the ensuing pressure gradient is able to counteract the gravitational acceleration, and a 
non-singular equilibrium state is reached. 

Mathematically, the evolution of a fluid with mass density ^(r, t) and velocity v(r, t) 
is described by the Euler equations, that amount to the conservation of mass 

dtp + V • Gov) = (1) 

and momentum 

VP 

d t v + (v • V) v = VU (2) 

P 

supplemented by the Poisson equation 

AU = 4ttG( P -(p)) (3) 

that relates the gravitational potential U(r, t) to the density distribution of the fluid, as 
well as by an equation of state that specifies the function P{p) or, equivalently, the sound 
speed 

*-% 

where P denotes the gas pressure. Finally, the initial data in p and v at some arbitrary 
time to have to be specified in order to close the system. 

The term (p) in equation (3) corresponds to the average gas density, and it vanishes 
for any system with finite mass contained in an infinite space. Here, we will consider the 
collapse of a protostar out of a gas cloud that is assumed to be roughly uniform on scales 
much larger than those relevant to the problem. Such conditions may be realized by a 
periodic cubic box of sufficiently large size L, i.e. a three-dimensional torus T 3 , as the 
spatial domain. In general, for a d-dimensional torus with periodicity L, 

(p) = Jdf Jd P&) dx - ( 5 ) 

We will also make the approximation, akin to Darcy's Law in Earth sciences, that 
the acceleration term that appears in the momentum conservation equation is always 
proportional to the fluid velocity 

d t v + (v • V) v « - (6) 
r 

where r is a constant with dimensions of time. Substituting in expression (2) and taking 
the divergence, one obtains 

-V • (pv) = -AP - V (pVU) , 

T 
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which, using the continuity and Poisson equations and writing the equation of state in 
terms of the sound speed, transforms into 

^ = c 2 s Ap + VpVU + AnGp {p - (p)) . 

T 

Since U is fully determined by p, our model consists on a single, non-linear and non-local, 
partial differential equation that describes the evolution of the gas density. 

In order to set the constant r, let us consider the collapse of a spherically-symmetric 
perturbation where the density decreases monotonically with radius. This will always be 
the case in the vicinity of a density peak, and, for Gaussian random fluctuations, it will also 
hold, on average, at large radii [!]. Under those conditions, the different spherical shells 
will not cross during the collapse phase, and the enclosed mass M(r) will be conserved for 
every shell. Starting from an initially homogeneous configuration, 



2GM 



2r 3 



|3tv + (vV)v| S^f V GM V 2tt<V 

so we adopt 

in order to recover the exact solution at arbitrarily small times. 

2 

Denoting the solution of Au = f — (/) by u = T(f), defining j3 = G J^ a L i , and choosing 
the units of time, length, and mass such that , 1 = 1, L = 2-ir, and (p) = 1, 

' b a /24ttG(p> 

respectively, we finally arrive at 

d t p = PAp + p(p-l)+Vp- VT(p) . (8) 

The case f3 = would be a model for a perfectly collisionless fluid (e.g. cold dark 
matter), whereas f3 > corresponds to an isothermal ideal gas. Moreover, we consider the 
generalized equation 

dtp = -f3A a p + p(p - 1) + Vp • VT(p) (9) 
where the fractional Laplacian A a = (— A) a / 2 is defined using Fourier theory: 

A^u= \£\ a u. 

Phisically, this is a nonlocal adhesion model. Due to the nonlocal character of the diffusive 
operators A a , the pressure forces now contain information about the total distribution of 
mass; in other words, the pressure at a given point x does not depend only on p(x) but 
on the overall distribution p(y) for all y. Such an equation may develop singularities in 
a finite time depending on the initial data. The case a = 1 is critical in the sense that 
we have a conservative derivative (the term Vp • Vi7) and a 'dissipative' one in the term 

Equation (8) can be re-written as 

d tP = f3Ap + V-(VUp), AU = p-(p), 

that bears obvious resemblances to the vorticity equation for the incompressible Navier- 
Stokes equations in two dimensions 

d t u) = vAuj + V • (V" 1 ^), A^ = w, 
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and it has been proposed by Jager and Luckhaus [20] as a simplified version of the classical 
(parabolic-elliptic) Patlak-Keller-Segel model [27, 22] of chemotaxis in biological systems 
(see e.g. [8, 10, 11, 12] and references therein; a recent review of the results concerning 
this equation can be found in [9]). The non-periodic (bounded or otherwise) case with 
(p) = has been previously studied in the context of gravitational collapse by P. Biler and 
co-workers, who have considered the problems of existence, conditions for blow-up, radial 
solutions, stationary solutions, and other qualitative properties [4, 5, 6]. L. Corrias, 
B. Perthame and H. Zaag have proved that, for small data in [|po||l<V2, there are global 
in time weak solution to equation (8) in d > 1, and blow up if the smallness condition 
does not hold (see Theorem 1.1 and 1.2 in [ ]). The periodic case has been studied by 
T.Senba and T.Suzuki in [ ]. In this paper they consider the case of measure valued weak 
solutions. 

We argue that there are several reasons why this equation constitutes a valid approach 
to the problem of gravitational collapse seeded by small perturbations of an initially ho- 
mogeneous medium in equilibrium (a situtation often encountered in Astrophysics). First, 
the family of hydrostatic equilibrium solutions of the full problem, 



are also equilibrium solutions of our approximation, since, taking the divergence of the 
above expression, 



In addition, the model is also consistent with more classical alternatives based on linear 
perturbation theory, being, in fact, very similar in spirit (though not exactly equivalent) to 
the Zel'dovich approximation [30] used in Cosmology to study the formation of large-scale 
structure in the primordial universe. In this approximation, particles move along straight 
trajectories of the form 



where r denotes the actual (Eulerian) position of the particle, q is its Lagrangian co- 
ordinate, the factor a(t) accounts for cosmic expansion, and the vector function p(q) is 
determined by the initial conditions. The evolution of the inhomogeneities is given by the 
growth factor b(t), chosen to match the results of Eulerian linear theory. For an Einstein- 
deSitter cosmology, which always provides a valid approximation for a matter-dominated 
universe at early times, b = a 2 , (^) = ^P-(p), and ^ = ^p(p)- Working in comoving 
coordinates, x = r/a, one obtains 



analogously to equation (7), although in this case {p) depends on time due to the expansion 
of space (see [14] for a detailed discussion of the cosmological implications of the parallelism 
between peculiar velocities and accelerations). 

The non-local generalization (9) has also received considerable attention during the 
past years (e.g. [19, 7, 13]). In particular, D. Li, J. Rodrigo and X. Zhang [23] have 
established local existence (using the linear semigroup and Duhamel's formula), as well as 
a continuation criterion derived by contradiction (see our Proposition 1) and some results 
concerning blow-up in finite time (Proposition 3). Global existence when the initial data 
is small in L 1 has been recently addressed in [25] and [26]. We will show that, in such a 
case, ||p||i/»(i) may be bounded uniformly, and global existence is obtained as a corollary. 
We state these results as a condition on the size of the pressure forces, /3; our sufficient 



VP + pVU = 



P Ap + p(p - 1) + Vp • VT(p) = 0. 



r(t) = a(t)q + 6(t)p(q) 
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conditions for gravitational stability are related to (and arguably supersede) the classical 
Jeans criterion, j3 > 1 [21], widely applied by the astrophysical comunity. 

The paper is organized as follows: In Section 2, we show that the problem is locally 
well posed (forward in time) in Sobolev spaces (Theorem 1) and a Beale-Kato-Majda-type 
criterion involving ||p||l°° (Proposition 1). These results have been proved previously (see 
Theorems 1.1 and 1.4 in [23]) by different techniques. We also present a result concerning 
the smoothing effect of our model (Theorem 2), due to the linear part when f3 > 0, and 
prove that the problem is ill posed backward in time. The smoothing effect result, up 
to our knowledge, is new but it is something that could be expected. The ill-posedness 
result is almost a consequence of the instant analyticity property. The conditions for 
gravitational stability, i.e. the decay of small perturbations (Theorem 3) and fluctuations 
of arbitrary size (Theorem 4), are investigated in Section 3. These results are compared 
with the well-known critical mass phenomenon in higher dimensions. Using Proposition 

1 these theorems implies global in time existence results which are, up to our knowledge, 
new for equation (9) (compare with Theorem 1.1 in [ ]). A finite time blow up result for 
some initial data in a special class is proved in Theorem 1.10 of [23]. These initial data 
are concentrated near the origin in terms of a parameter 5 such that, for 5 « 1, the 
maximum of the initial datum grows, roughly speaking, like 1/5. Taking this Theorem 
into account, our global in time result for small data seems to be quite sharp. In this 
Section we also give a minimum time for blow up which is the periodic version of Theorem 
1.9 in [23] . Numerical simulations are presented in Section 4. 

2 Well-posedness 

2.1 Well-posedness in Sobolev spaces 

For the sake of mathematical consistency, we show in this section that equation (9) is well 
posed if /3 > and < a < 2. First, let us note that the average density (p) defined in 
equation (5) is conserved during the evolution of the system, and thus it is a constant that 
only depends on the initial data 1 . We state this claim in the following Lemma: 

Lemma 1 (Mass conservation). Let p(x,t) be a classical solution. Then, the total mass 
is conserved 

/ p(x,t) dx = I p(x,0) dx, 

Jjd. Jf d 

and, as a consequence, (p) is a constant. 

Proof. We observe that, using 2-7r/(0) = J T f(x) dx, the term 

/ A a u = 0. 

Now, integrating and using the Divergence Theorem the proof follows. □ 

To prove the existence and uniqueness of classical solution, we proceed as in [3]. We 
obtain some 'a priori' (i.e. assuming the existence of classical solution) bounds for the 
usual norm in the space H k (T d ) and then regularize equation (9) so that all the regularized 
problems have a classical solution for a uniform time T. To conclude, we use the 'a priori' 
bound to show that the solutions to the regularized problem form a Cauchy sequence 
whose limit is the solution to the original equation. In order to simplify the notation, we 
will abbreviate p(x,t) = p(x), or simply p, throughout the rest of the paper. 

1 Without loss of generality, one can set (p) = 1 by an adequate choice of units. 
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Theorem 1 (Well-posedness) . Given equation (9) with (3 > 0, < a < 2, and initial 
data < po(x,0) G H k (T d ), where k > d/2 + 2, there exists a time r = r(po) > 
and a unique solution p £ C([0,t],H 8 ) n L°°([0, r], H k (T)), for s < k. Moreover, p £ 
C 1 ([0,r] 1 C(T))nC([0,r],C 2 (T)). 

Proof. For simplicity, we only work out in detail the case k = 3 in one spatial dimension, 
the general case with d = 2, 3 and k > 3 being analogous. 
(Existence:) We define the norm in H k (T) as 

II II 2 — II II 2 _l II ,a* II 2 
II " \\H k — II ' Hi 2 + ll°-r ' Hl 2 - 

We multiply the equation by p and integrate by parts: 

p{x)d tP {x)dx = -p f p{x)A a p{x)+ [ (p(x)) 2 (p(x) - l)dx+ [ 9x( " P [ X ^ 2 d x Udx 



Using Parseval Theorem and Holder inequality we get the estimate 

j t \\ P \\ 2 L 2 < -2p\\A a/2 p\\h + ci\\ph°°\\p\\h < <\\p\\h» + i) 2 IIpIIh3. (10) 

In the same way, for ||^p||? 2 , we obtain the following bound: 

^ll^llia < c(\\p\\L*> + ||^IU-)l|p|llfa < 4p\\h-4\\p\\h3 + I)'- (11) 
Using (10) and (11), one finally arrives to 

^.\\p\\H-i < c(\\p\\ H 3 + I) 2 , 

so 

ii || ^ WpoWh 3 / 10 n 
WpWhs < i n — n — t- (12) 

We only sketch the rest of the proof because it is classical. In order to show the 
existence of classical solutions to problem (9), we consider J £ C^°, J(x) = J7"(|a;|), 
J > 0, and L J = 1. For e > 0, we define 

J e (x) = I J (^) (13) 
and consider the regularized problems 

dtp" = -PA * ^Je ^P' + Je* {J, * P%Je * P* ~ 1)) + Je * (V ' & * p € ■ VT(^)), 

where 

AT(p e ) = (J- e *p e -l). 

In these regularized problems, the total mass is also conserved, and the previous bounds 
(10) and (11) hold. Using classical energy methods (see [3]), one may prove that the 
regularized problems have an unique smooth solution and, using expression (12), that the 
sequence p e , indexed in e, is Cauchy in C([0, r], H s ), < s < 3, with r = r(po) < cUpoTT' 
Thus, it has a limit p, which is a classical solution of our problem. In order to show that 
p € H 3 we use that p e is uniformly bounded in H s and, due to the strong convergence to 
p in H s with s < 3, this implies that p e are weakly convergent to p in H 3 . 
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(Uniqueness:) Suppose that for the same initial data po we have two classical solutions 
Pi, P2 £ C([0, t], H 3 (T d )). Then the equation for p = p2 — p\ is 

d t p = -(3A a p+( P2 (p 2 - 1) - pi(pi - 1)) + Vp 2 • VT(p 2 ) - Vpi • VT(pi), 

so 

~-j>ll£s < -£l|A a/2 p||£ 2 + \\p\\U\\P2\\l~ + \\pi\\l°°)+ / pV Pl V(r(p 2 )-r( Pl ))dx. 
2 est J T d 

We observe that the equation for T(p 2 ) — T(pi) is 

A(T(p 2 ) - T( Pl )) = p. 
Using Poincare inequality we have that 

\\V{T( P2 ) - T{ Pl ))\\ L 2 <c\\p\\ L 2. 

Then 

/ pV Pl V(T(p 2 )-T( Pl ))dx< \\Vp 1 \\ L ~\\p\\ L 2\\V(T(p 2 )-T(p 1 ))\\ L 2 <c||Vpi||^||p||| 2 . 
Putting all together we obtain 

j t \\p\\h ^ c (P2,pi)||p|||2- 
Applying Gronwall inequality we conclude the uniqueness. □ 

We give an alternative proof of a Beale-Kato-Majda-type criterion characterising the 
possible singularities (see [23]): 

Proposition 1 (Singularity formation). Suppose that we have 

•t+ 

\\p\\l°° < oo, 

Then the solution exists upon time t + + a for a sufficiently small a . 

Proof. To clarify the exposition, we consider the case with d = 1 first. To prove the result 
we will assume that + 

/ \\p\\ L oa(s)ds = M < oo (14) 
J o 



Using (10) and (11) we have 

jtJHIh 3 ^ c (\\p\\l°° + II^p||l°o + l)||p||^3, 
so, using Gronwall inequality we conclude that 

\\p\\m(t) < \\p \\ H3 e c ^\\p\\^^)+\mp\\ L ^(s)ds+t_ 

Now we use the technique exposed in [17]: we use Rademacher Theorem for d 2 p(xt) = 
W^xPWl 00 (£)■ We can do that if the solution is sufficiently smooth. The evolution of this 
quantity is given by 

d t d 2 x p(x t ) = -pA a d 2 x p(x t ) + d 2 x p(x t )(3(p(x t ) - 1) + p(x t )) + 3(d xP (x t )) 2 . 
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We use the generalized Landau inequality (see [24]) 

\d x p(x)\ 2 <2p(x)\\d 2 x p\\ L ~. 

Due to this pointwise inequality we observe that the quadratic terms (d x p(xt)) 2 are linear 
in d 2 p(xt). Thus we obtain 

d t d 2 p(x t ) < W\\p\\ Lao (t)d 2 x p(x t ) ||^|U-(t) < ||^po|Ucce 10 J"o 10 IWI^W*-, 

and we conclude the result in the one dimensional case . 

Now we consider the case d = 2. The first step is to obtain a bound on ||Ap||^2. We 
multiply the equation by A 2 p and integrate in space. The dissipative terms are negative 
so we neglect it. The reaction term contributes with 

/ Ap(x)A(p(x)(p(x) - l))dx < c\\Ap\\l 2 \\p\\ L oo + ||Ap|| l2 |||VH 2 ||l 2 , 

and using a classical Gagliardo-Niremberg interpolation inequality we get 

|||Vp| 2 || i2 < c\\Vp\\ 2 L , < c\\Ap\\ L2 \\p-l\\ L oo. 
Thus the contribution of the reaction term can be bounded as 

/ Ap(x)A{p(x){p(x)-\))dx<c\\Ap\\ 2 L 4p\\ L ~. (15) 

The contribution of the transport term is 

/ A 2 p(x)Vp(x) ■ VU(x)dx = Jx + J 2 + J 3 + J a + J5, 

where 

Ji = / Ap{x)VU(x) ■ VAp{x)dx < c||A/o[|?2[|/o[|l">, 
Jt 2 

h = [ Ap(x)\Vp(x)\ 2 dx < c\\Ap\\ L 2\\\Vp\ 2 \\ L 2 < c\\Ap\\l 2 \\p\\ L ~, 
Jt 2 

J 3 = 4 f Ap{x)d 2 p{x)d 2 U {x)dx < c||Ap|| 2 2 1|0 2 , C/|| L ~, 
Jt 2 

J 4 = 2 / (Ap(x)) 2 d 2 Xl U(x)dx < c\\Ap\\ 2 L2 \\d 2 xl U\\ L ~, 

J 5 = 2 I Ap{x)d 2 X2 p(x){p{x) - l)dx < c||Ap|| 2 2 ||p|| L ~, 
Jt 2 

and, putting all together we obtain 

j t \\M\L*< c\\Ap\\ L2 (\\p\\ L ~ + \\d 2 Xl U\\ L ~). 
Due to Gronwall inequality we obtain 

\\p\\ H 2{t) < c{M)\\p \\ H 2 exp H Wd^UW^is) + \\d Xl d X2 U\\L°°(s)ds^ . (16) 

We use a classical inequality (see 3.2b and Appendix 1 in [28] and [2]) 

\\V\\ L °° < c(||C||lp + IICIU- max{l, log(||C|| w .-WHCl|L»))}, l<p<oo,s>l + 2/p 
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where V(C) is a singular integral operator of Calderon-Zygmund type of £. In our case we 
have 

\\d Xl d X] U\\ L oo < c(\\p - 1\\ L 2 + \\p - l\\ L oo max{l,log(||p - l||#2/ll/° - 

< c(M, ||po||L2)(l + \\p\\l°° max{l,log(c||p|| H 2)}). (17) 

Inserting in (16) we have 

IMIWt) < c(M)\\p \\ H 2 exp (c(M, Upollxa) J (1 + ||p||i/x,(s)log(c||p||#2(s) + e)ds) . 

Now let y{t) = c\\p\\ H 2(t) + e. We have 

y(t) < c(M)y(0)exp (c(M,\\ Po \\ L 2) (1 + ||p||^(s)log(y(s))ds^ , 

and write z(t) = log(y(t)). We get 

z(t) <c(M) + z(0) + c{M,\\p \\ L 2)t+ [ \\p\\ L oo(s)z(s)ds, 

Jo 

and due to the integral version of the Gronwall inequality we get 

z{t) <c{M,\\p \\ H 2)(l + t). (18) 

Thus, we conclude that the bound (14) implies p(t) G H 2 . To obtain higher regularity 
we can do the same. Indeed, taking derivatives and multiplying by the correct terms we 
obtain 

\\p\\h* < WpoWh 4 ex P {^j + II^Plli 2 ( s ) + ll^xi^lli 00 ^) + \\dx 1 X2 U\\ L ^(s) S j ds, 

where we used the following classical interpolation inequalities 

\\Dff L4 <c\\D 2 f\\ L 4f\\ L ~, and \\D 2 f\\ 2 L , < \\D 2 f\\ L2 \\D 2 f\\ L ~. 
Using the previous bounds (14) (17) and (18) we conclude the result. 



□ 



2.2 Smoothing effect 



Our proof follows the steps presented in [16]. The proof is based in a priori estimates for 
the complex extension of the function p on the complex strip S = {x + it;,x G T, |£| < kt}, 
for certain constant k > (see Figure 1). We define 



IMIl2(S) = J^\p{x + ikt)\ dx + J \p(x-ikt)\ dx, \\p\\ H 3 {s) = \\p\\ L 2 (s) + \\d x p\\ L2{s) . 

For the complex extension the equation is 

d t p(x + i£) = -f3A a p(x + i£) + p(x + i$)(p(x + i£) - 1) + Vp(x + if) • VT(p)(x + if) (19) 

Theorem 2 (Smoothing effect). Let us consider equation (9) with a > 1 urti/i po G -ff 3 
as initial data and T as spatial domain. Then the classical solution p (which exists, at 
least locally in time, because of Theorem 1 ) continues analitically into the strip S for times 
0<t<r(p ). 
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t t 




Figure 1: The strip of analyticity for p. 

Proof. For the sake of brevity we work with both terms ±ikt at the same time. Firstly 
we consider the critical case a = 1. We start with the L 2 (S) norm: 

— — I \p(x ± ikt)\ 2 dx = Re f p(x ± ikt)(dtp(x ± ikt) ± ikd x p(x ± ikt))dx. 
2 at J j Jj 

Using Plancherel Theorem we have 

-P / p(x ± ikt)Ap(x ± ikt)dx = -(3 / \A 1/2 p(x ± ikt)\ 2 dx < 0. 
Jt Jt 

The reaction term can be bounded as follows: 

Re J p(x ± ikt)p(x ± ikt)(p(x ± ikt) — l)dx < (||p||l°°(s) + l)ll/ 9 lli 2 (5)- 



We have 



I\ = Re J p(x±ikt)d x p(x±ikt)d x U(x±ikt)dx = Re / <9zl^(x±ifci)^<9 x |/9(x±i£;i)| 2 dx 

+ Re i / c^t/ (x ± ikt)Rep(x ± ifci)<? x Imp(a; ± ikt)dx 
Jt 

— Re i d x ll (x ± i/ct)Imp(x ± ikt)d x Rep(x ± ikt)dx. 
Jt 

Integrating by parts we have 
Ji = Re A d x ll(x ± z/ci)-<9 x |p(x ± ikt)\ 2 dx 

JT 

+ 2 / c^imt/^x ± i/ct)Im / o(x ± ikt)d x Rep(x ± ikt)dx, 
Jt 
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and we have 

h < 4p\\h(s)(\\p\\L°°(S) + 1) + ||3 a! Iml7|| L o (s)||p|| £ 2( S )||a a .p|| £( 2( S ). 

Due to the mean conservation and that the initial data has Impo = we have the bound 

\\d x ImU\\ L oo( S ) < cll^ImJ/II^^II^ImJ/H^ 2 ^ < c||Imp|| L 2(s), 

where in the last inequality we use the Poincare inequality. 
The last term can be bounded as follows: 



Re / p(x±ikt)(±ikd x p(x±ikt))dx < k\\p\\L2(s)\\dxP\\L 2 

J 7 

Putting all together and using Sobolev embedding we have 

^||p|lia( S ) < c||p||^3 (5) (l + ||p||if3 (5) ). 



(5)- 



(20) 



We have 

ld_ 

2di JT 

Following the previous techniques we have 



\d x p(x±ikt)\ 2 dx = Re / d z x p(x ± ikt){d t dlp{x ± ikt) ± ikd\p{x ± ikt))dx. 

Jt 



1 d 
2~dt 



[ \dlp(x±ikt)\dx<c\\p\\% 3{s) (\\p\\ H s {s) + l)-P [ \k l/2 dlp\ 2 dx 

JT JT 

/ d x pd x pd x lldx ± kRe i / d x pd x pdx. 
Jt Jt 



+ Re 



We write 



J\ = Re / d x pd x pd x lldx and J2 = ±kRe i / d x pd x pdx. 
Jt Jt 

Splitting into real and imaginary part we have 

Ji = Re f d x u\d x \dlp\ 2 dx+ [ d^Repd x lmplmd x Udx - [ d^lmpd x Replmd x Udx 
Jt 2 Jj Jj 

<c||^p|| L 2 (5) (||p|| L oc (5) +l) + 2K 1 , 

where the last integral is 



dlRepAHd x lmplmd x Udx 



A 1/2 HdllmpA^ 2 (dlReplmd x U)da 

< IIA^^pll^^HA^^RepIm^i/)!!^^. 

We use the commutator estimate 

\\A 1/2 (fg) - fA 1/2 g\\ms) < c||^/|| L oc (5) || 5 || L2(5) 

to conclude that 

K x < HAV^^pll^^dl^ImC/IUoo^ll^pll^^ + II^ImC/IUoo^HA 1 / 2 ^!!^^)) 

< c(||AV 2 ^ p || La( ^||p|| Loo(s) ||^ p || La( ^ + IIImpll^^HA 1 / 2 ^!!^^). 
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If we use Young inequality we have 

K x < (| + c||Imp|| i2(s) ) || k x / 2 ^ xP t L2{s) + ^|Mlioo (s) ||^p||£ a(s) . 
We can also bound 

J 2 < k\\A l / 2 dl P \\ L 2 {s) . 

Finally, taking e = /3/2 
d f 

— \dlp{x ± ikt)\dx < cllpll^s^) (||p||h3 (5) + 1) 

+ (c||Imp|| L2( 5) +k-^fj \\A l / 2 ^p\\ L2{S) , (21) 

and this is a correct bound if the term 

3/3 

K 2 (t) = c\\lmp\\ L 2 {s) +k--^<0. 

Initially 1^(0) = k — so taking k = /3/4 < /3 the condition K 2 < will be satisfied for 
a (maybe short) time if we can bound the evolution of ||Imp||£2(5). 
We define the 'energy' 



\p\\a(s) = Wp\\h3(s) + — 



2 c|l Im PllL2(S) 

The evolution of this energy is 



^IIpIU(S) < c ((IIpIU(S) + !) 2 + \\p\\A(s)^l\\ lm P\\L^S))- 
At this point it is easy to obtain 



^\\ lm P\\L^S) < 4\\p\\h3(S) + 1 )ll Im / 5 lli2(S)' 

so 

J t \\p\U(S) < (IIPIU(5)+1) 5 

and we conclude that there is a time r{p) so that ||/o||a(s) < 00 • 

Now we approximate the problem in the real line as in the proof of the Theorem 1, 
using the rescaled heat kernel as a mollifier: 

dtp' = -0 'Je * A a Je *p e + Je* (Je * p\Je * p" ~ 1)) + Je * (VJ e * p* • VJ e * T(p £ )), 

where 

AT(p e ) = (J e *p e -l), 

and we use p$ = J e * po as initial data. Using Picard's Theorem in H S (T) we obtain 
solutions p e up to time r € which are analytic. Due to the above estimates we have a time 
a common time of existence r depending on the initial data po but not in e and if t < r(po), 
\\p £ \\h 3 (S) < 00 • Now we pass to the limit in the Hardy-Sobolev space H 3 (S) obtaining 
p and we use the fact that if the energy || • \\a(s) remains bounded then this implies the 
analyticity. Due to the uniqueness, this is the same p as in Theorem 1. 
In the proof for the case a > 1 we use the inequality 

l|A 1/2 T||i 2( 5) < ||r||| 2( 5 } + ||A"/ 2 T||i 2(5) 

obtained by splitting in the Fourier side. We conclude the proof in the supercritical case 
following the same steps. □ 
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Using the smoothing effect, one can prove ill-posedness in Sobolev spaces when the 
problem is considered backward in time. We remark that this result does not require 
global existence of solutions. 

Proposition 2 (Ill-posedness). There are solutions p to the backward in time equation 
(9) with /3 > 0, such that ||/5||#fc(0) < e and \\p\\fj k (^) = 00 f or a ^ < e, sufficiently small 
< 5, and k > 3. 

Proof. Let us fix k = 4, the other cases being similar. Take go(x) £ H^(T) but go ^ H 4 (T) 
and consider the solution (forward in time) p x to the equation (9) with initial data p(x, 0) = 
Xgo(x) where < A < 1. We have (see Theorem 2) that there is a uniform (in A) time 5* (go) 
such that p x exists and is analytic up to time 5* (go)- Now define p X ' S (x, t) = p x (x, —t + 5) 
for fixed < 5 < 5* (go). We have that Hp^H^^) = A||go||.H" 4 = oo. Now, using Theorem 
2 we have that there is a growing strip of analitycity for p x where we can apply Cauchy's 
integral formula in order to obtain 

ll^ 4 ^Vi| (n\ - \\Mn x \\ (X\ <r C ^ iifl3„An / c(f3,g )X 
WOxP IIl 2 (t)(U) - \\o x p \\ L 2(t)(o) < \\ d xP Wl 2 (s s *) < ^ • 

Taking < A < min{l, e6/c(/3,go)} we conclude the proof. □ 



3 Gravitational stability 

The blow-up of the solutions in finite time is of particular interest from the physical 
point of view, since it determines whether a gravitationally bound structure may form 
or not. In order to investigate the necessary conditions for star formation, the present 
section is devoted to the evolution of ||p||z,°° (*)- More precisely, it is our goal to determine 
whether /when the growth of singularities can be avoided in our non-local adhesion model. 
We will focus on the critical case (a = 1) and remark that the resulting equation is more 
dissipative if a > 1. All our results are expressed in units such as (p) = 1; the conversion 
to other systems can be trivially carried out by substituting every occurrence of p by -^y 
or, alternatively, (3 by fi{p). 



3.1 Minimum time for star formation 

First of all, let us address the minimum time required for the gravitational collapse of a 
protostar. For those cases where a star is formed, i.e. p(x*,t*) — > oo, we provide a lower 
bound for the value of t* as a function of the initial condition. Note that, by virtue of 
Proposition 1, classical solution is guaranteed up to a time £*. 

Proposition 3 (Minimum time for blow-up). Let p > be a smooth solution of (9) in 
T d , d < 3. Then we have 



dt 
thus 



1 "p\\L°°(t) < ||p||z,°°(i)(|M|i,<*>(t) - 1), 



i ii n\ ^ Po i° 



Ipo||l°° + (1 - ||po|U°°)e*' 
As a consequence the minimum time for the formation of a star is 



t* = log 



IpoI|l°° — 



13 



Proof. Due to the smoothness of p and using Rademacher Theorem, we have that the 
evolution of 

= raaxp(x,t) = p(x t ), 

X 

is given by 

j t p{x t ) < -PA a p(x t ) + p{xt){p{x t ) - 1) < P(xt)(p(x t ) - 1), 

where the last inequality is due to the integral representation for the diffusive operator. 
Now the result is straightforward. □ 

3.2 Exponential suppresion of density fluctuations 

Our main result, and the one that bears more physical relevance, is that equation (9) with 
a = 1 does not have a blow-up for ||p||l°° (t) - i.e. no star can ever form - if (3 is big enough. 
Physically, f3 is the dimensionless ratio of the free-fall time, representative of the strength 
of the gravitational forces, and the sound-crossing time, representative of the restoring 
pressure force. Depending on the amplitude of the initial conditions (more precisely, the 
value of ||/Oo||l°°), there is a critical value of the sound speed above which the pressure 
gradient is able to overcome gravity. Or, in other words, for a uniform gas cloud at a 
given temperature (i.e. /?), sufficiently small perturbations cannot grow. Mathematically, 
this result is a Maximum Principle for ||p||l°° for a special class of initial data, and the 
description in terms of f3 or HpoIIl 1 merely reflects the choice of units. Both choices are, 
of course, absolutely equivalent. 

Consistently with the results of numerical simulations, we find that, for sufficiently 
small perturbations and/or large pressure forces, there is a depletion of the solution (see 
Figure 2), and the system evolves towards the homogeneous equilibrium state p(x) = 1. 
Since ||p||i°° is bounded, global existence of classical solution follows from Proposition 1. 




-3 -2 -1 1 2 3 -1 



Figure 2: a) The initial data (in blue) and the evolution (in red) for the one-dimensional 
case of equation (9) with /3 = 1 and b) the complete evolution. 

Theorem 3 (Decay of small perturbations). For (9) with < a < 2, spatial domain T d , 
d = 1,2, 3, and initial data po, if 

\\po\\l°° < c a ,d(3, 

where c a ^ is an explicit constant, then the maximum density decays exponentially, 

\\ p \\ LOO (t) < i + [ UpoIU— — i] e-K^-i^M*. 
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Figure 3: a) The initial data (in blue) and the evolution (in red) for the one-dimensional 
case of equation (9) with /3 = and b) the total evolution. 



Moreover, we have that ci.i = 1. 

Proof. First, we prove the case a = 1, d = 1. Denoting the point where p achieves its 
maximum value by xt and using the integral representation of the fractional Laplacian, 
one can easily derive the bound 

A„(z,) - ip.v. jf l^jjffi, > I / T pM - ,M iy - ?M - i 

Thus 

d t p(x t ) < -P[p(xt)-l]+p(xt)[p(x t )-l] = [p(x t )-l][p(x t )-P], 

and the maximum density will always decrease if p(xt) < /3. In fact, 

d t [p(x t )-l] < - [/3-\\p \\ L ~}[p(x t )-l]. 

Now we proceed with the d— dimensional case. From the integral representation of the 
fractional Laplacian in T d 



A a p(x t ) = C aA V P.V. / 



p{xt) - p(y) , 
\x t -y- 2irv\ d + a 



we have that 

Ay>C ad F.V. [ "M-^ dy. 
' Jw \xt-y\ d+a 

Since \x t — y\ < nyd, 

A°P > c a4 [p(x t ) - 1} . 

and we can conclude the result analogously to the one dimensional case. □ 

In one spatial dimension, it is also possible to prove that, for sufficiently large f3 
(or, equivalently, sufficiently small || poll 2,1)1 a ^ density fluctuations are exponentially sup- 
pressed, regardless of their initial amplitude. 

Theorem 4 (Global stability for d = 1). Given equation (9) with a = 1, spatial domain 
T, and initial data po, if 

P > 4vr 2 , 
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then the maximum density decays exponentially, 

\\p\\L°°{t) < 1 + [ ||po||z°° - 1] e~*. 
Moreover, if 1 < a < 2 then, for any f3 > 0, the following inequality holds 

\\p\\l°° < c(a,(3). 

Proof. Let us consider the neighbourhood around the absolute density maximum 

= {y I \x t -y\<25 }, 

where 5 is a parameter that will be set later, and decompose it in two subsets, Q$ = U\\jlA<i-, 
such that 

Ui = { y G n 5 I p(x t ) - p(y) > p(x t )/2 }, 

and 

U 2 = { y G Oj I p(x t ) - p(y) < p(x t )/2}. 

Since |x| > | sin(x)|, 



p(s t ) - p(y) , . 1 
— ay > — 



Substituting | S~2^ | = 4<5 and using that 



Wi 



<5 2 



1 g(g«) 
2vr 2<5 2 



(|fi«|-|W2|). 



2vr 



one obtains the lower bound 

Ap(x t ) > 
which is maximized for the choice 



/ p(y)dy 
Jt 

1 p{xt) 



>> 17-/ I 

> — 1^2 1, 



2^ 25 2 
2tt 



45 



4"7T 



P(»t)' 



It is important to note, though, that the compactness of the domain also imposes that 
5 < 7r/2 or, in other words, 

p(x t ) > 4. 

Inserting the bound for the dissipative term in equation (9), we have 

P 2 ( x t) , \ r ^ 1 1 _2,_ x f 1 Z 3 



< -(3 ■ 



4tt 2 



+ p(x t ) [p(^t) - 1] = P (xt) 



1 



4vr 2 



p(xt) 



UP> in 2 , 

dtp(xt) < -p(x t ) < -[p(x t ) -1], 

which proves the desired result for p(x<) > 4. Theorem 3 implies that this bound (indeed, 
a stricter one) also applies for p(xt) < 4. This concludes with the first part of the result. 
In the case 1 < a < 2 we have 



Ap(x t ) > 



1 p{xt) 
2vr 25 1+a 



46 



47T 

P{xt) 



Consequently, we get 



<■>//'<■'■/):£ -P^^+p(x t )[p(x t )-l] 



and we conclude the proof. 



□ 
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Remark 1 Recall that in 2D (8) has a critical mass phenomenon, i.e. if ||po||l°°(r 2 ) < 8vr, 
then there is a global solution, while if ||po||l°°(r 2 ) > 8ir there is a finite time blow up. If 
we take /3 = 1 and an initial datum with arbitrary L 1 norm, Theorem 4 gives us that, if 
H/00 Hl 1 < 2^ then the maximum decays exponentially and the solution is global. 

4 Numerical simulations 

In order to simulate the equation (9) we use a Fourier-collocation method (see [ 1 5] ) . We 
use Fast Fourier Transform techniques to solve the Poisson equation and to approximate 
the spatial part of (9). For the time integrator we use an explicit Runge-Kutta of order 4. 

We center our atention in the 2D case. The other cases are analogous. We consider an 
equispaced discretization of our spatial domain T 2 , 

TN = {(x k ,yi) ■ x k = ^(2k-N), yi = ^(2l-l),k,l = 0,...,N-l}, 

where N is the number of points on each interval [— 7r, it]. We consider as our approximate 
solution the pair p(x,y) and T(p)(x,y) in the grid We have 

N/2-l 

p(x k ,yi)= Pn, m e inXk ^ mm , -{n 2 + m 2 )T{p) nm = ^Gp- n ^ T(p) 00 = 0, 

n,m=-N/2 

and it is important to remark that the coefficients p n ^ m are related to the so-known Fourier 
coefficients, typically p n ,mi brtt. they are not the same (see [15]). In the same way we can 
approximate all terms in (8). Once we compute derivatives of p and U via the FFT in the 
equation, we get back to the physical space using the IFFT. 

This ends the space discretization part. We consider our equation in the time interval 
[0, t], for some r > 0. As before, let Qj an equispaced partition of the time interval, 

Qj = {t s : t s = st I J, s = 0, . . . , J - 1}, 

where J is the number points in the time grid. For the evolution part of the equation, 
dtp, we are going to use a Runge-Kutta method (specifically, the classic explicit method 
of order 4), and we start with the initial data 

po(x, y) = | sin x\ + | sin y \ 

represented in Figure 4. 

We studied the evolution of IHIl 00 ^) for various f3 G [0,1]. On the left panel in 
Figure 5, we highlight in red the case /? = 0.2 and in blue f3 = 0.8. On the right panel, we 
plot the results for values of f3 in the interval [0.32, 0.34]. 

The mechanism for singularity formation, i.e. the blow-up for HpIIl 00 ^)) is illutrated 
in Figure 6. In Figure 7, one can see how the solution is depleted for large (3. 
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